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Abstract 

We study instanton contributions to hadronic decays of the scalar glueball, the pseudoscalar 
charmonium state rj c , and the scalar charmonium state Xc- Hadronic decays of the rj c are of par- 
ticular interest. The three main decay channels are KKir, nirn and n'inr, each with an unusually 
large branching ratio ~ 5%. On the quark level, all three decays correspond to an instanton 
type vertex (cc)(ss)(dd)(uu). We show that the total decay rate into three pseudoscalar mesons 
can be reproduced using an instanton size distribution consistent with phenomenology and lat- 
tice results. Instantons correctly reproduce the ratio B(inrn) / B(irirr]') but over-predict the ratio 
B(KKtt)/ B(inrn(rj')). We consider the role of scalar resonances and suggest that the decay mech- 
anism can be studied by measuring the angular distribution of decay products. 
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I. INTRODUCTION 



The charmonium system has played an important role in shaping our knowledge of per- 
turbative and non-perturbative QCD. The discovery of the J/ip as a narrow resonance in 
e + e~ annihilation confirmed the existence of a new quantum number, charm. The analysis 
of charmonium decays in e + e~ pairs, photons and hadrons established the hypothesis that 
the J/ip and r\ c are, to a good approximation, non-relativistic 3 Si and 1 S bound states 
of heavy charm and anti-charm quarks. However, non-perturbative dynamics does play an 
important role in the charmonium system [1, 2]. For example, an analysis of the ip spectrum 
lead to the first determination of the gluon condensate. 

The total width of charmonium is dominated by short distance physics and can be stud- 
ied in perturbative QCD [3]. The only non-perturbative input in these calculations is the 
wave function at the origin. A systematic framework for these calculations is provided by 
the non-relativistic QCD (NRQCD) factorization method [4]. NRQCD facilitates higher 
order calculations and relates the decays of states with different quantum numbers. QCD 
factorization can also be applied to transitions of the type ip' — > ip + X [5, 6]. 

The study of exclusive decays of charmonium into light hadrons is much more complicated 
and very little work in this direction has been done. Perturbative QCD implies some helicity 
selection rules, for example i] c -j^ pp,pp and J/ip pn, pa\ [7, 8], but these rules are strongly 
violated [9]. The J/vp decays mostly into an odd number of Goldstone bosons. The average 
multiplicity is ~ (5—7), which is consistent with the average multiplicity in e + e~ annihilation 
away from the J /ip peak. Many decay channels have been observed, but none of them stand 
out. Consequently, we would expect the r\ c to decay mostly into an even number of pions 
with similar multiplicity. However, the measured decay rates are not in accordance with 
this expectation. The three main decay channels of the r\ c are KKit, rjirir and rj'mr, each 
with an unusually large branching ratio of ~ 5%. Bjorken observed that these three decays 
correspond to a quark vertex of the form (cc)(ss)(dd)(uu) and suggested that r\ c decays are 
a "smoking gun" for instanton effects in heavy quark decays [10]. 

In this paper we shall try to follow up on this idea by performing a more quantitative 
estimate of the instanton contribution to rj c and Xc decays. The paper is organized as follows. 
In section II we introduce the instanton induced effective lagrangian. In the following sections 
we apply the effective lagrangian to the decays of the scalar glueball, eta charm, and chi 
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charm. We should note that this investigation should be seen as part of a larger efffort 
to identify "direct" instanton contributions in hadronic reactions, such as deep inelastic 
scattering, the AI = 1/2 rule, or r\ production in pp scattering [11-14]. 

II. EFFECTIVE LAGRANGIANS 

Instanton effects in hadronic physics have been studied extensively [15, 16]. Instantons 
play an important role in understanding the U(1)a anomaly and the mass of the rf . In 
addition to that, there is also evidence that instantons provide the mechanism for chiral 
symmetry breaking and play an important role in determining the structure of light hadrons. 
All of these phenomena are intimately related to the presence of chiral zero modes in the 
spectrum of the Dirac operator in the background field of an instanton. The situation in 
heavy quark systems is quite different. Fermionic zero modes are not important and the 
instanton contribution to the heavy quark potential is small [17]. 

This does not imply that instanton effects are not relevant. The non-perturbative gluon 
condensate plays an important role in the charmonium system [1, 2], and instantons con- 
tribute to the gluon condensate. In general, the charmonium system provides a laboratory 
for studying non-perturbative glue in QCD. The decay of a charmonium state below the 
DD threshold involves an intermediate gluonic state. Since the charmonium system is 
small, r cS ~ (fm c ) _1 < Aqq D , the gluonic system is also expected to be small. For this 
reason charmonium decays have long been used for glueball searches. 

Since charmonium decays produce a small gluonic system we expect that the cc system 
mainly couples to instantons of size p ~ r c5 ~ {vm c )~ l . In this limit the instanton effects 
can be summarized in terms of an effective lagrangian [18-21]. 



where t a = \\ a with tr[A a A b ] = 25 ab are SU(3) generators, 1 2 = diag(l,l,0), rf^ is the 



coordinates, the instanton position z, the instanton size p, and the color orientation U G 
SU(N C ). We also define the rotation matrix R ab by R aa ' X a ' = UX a Ul For an anti-instanton 




(1) 



't Hooft symbol and a 



ibV'T*']- The instanton is characterized by AN C collective 
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we have to replace L <-> R and rj <-> rj. The semi-classical instanton density d(p) is given by 

--^1 (2) 

where g(p) is the running coupling constant. For small p we have d(p) ~ p b ~ 5 where 
b = (lliV c )/3 — (2Nf)/3 is the first coefficient of the beta function. 

Expanding the effective lagrangian in powers of the external gluon field gives the leading 
instanton contribution to different physical matrix elements. If the instanton size is very 
small, p <C TnJ 1 , we can treat the charm quark mass as light and there is an effective 
vertex of the form (uu) (dd) (ss) (cc) which contributes to charmonium decays. Since the 
density of instantons grows as a large power of p the contribution from this regime is very 
small. In the realistic case p ~ (t>m c ) _1 we treat the charm quark as heavy and the charm 
contribution to the fermion determinant is absorbed in the instanton density d(p). The 
dominant contribution to charmonium decays then arises from expanding the gluonic part 
of the effective lagrangian to second order in the field strength tensor. This provides effective 
vertices of the form (GG) (uj^u) (dj^d) (s^^s) , (G 2 )(ifj 5 u)(d'j 5 d)(ss), etc. 

We observe that the Nf = 3 fermionic lagrangian combined with the gluonic term ex- 
panded to second order in the field strength involves an integral over the color orientation 
of the instanton which is of the form J dU^UijU^) 5 . This integral gives (5!) 2 terms. A more 
manageable result is obtained by using the vacuum dominance approximation. We assume 
that the coupling of the initial charmonium or glueball state to the instanton proceeds via 
a matrix element of the form (0 ++ |G 2 |0) or (0~ + \GG\0). In this case we can use 

<0 ++ |G^G» = 12{N \_ 1 f b (^^ - ^a){0 ++ \G a ; a G a ; a \0) (3) 

in order to simplify the color average. The vacuum dominance approximation implies that 
the color average of the fermionic and gluonic parts of the interaction can be performed 
independently. In the limit of massless quarks the instanton (I) and anti-instanton (A) 
lagrangian responsible for the decay of scalar and pseudoscalar charmonium decays is given 

by 

C I+A = jdz^dp {N f_ P l )as { {G 2 ~ GG) x L ftI + (G 2 + GG) x L ftA } . (4) 
Here, £f,iA is the color averaged Nf = 3 fermionic effective lagrangian [15, 16, 21]. 
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p 5 (jV c -l)!(JV c -2)! _ W 2 



„ x d (p) 0.466 exp(-1.679iV c )1.34 7V / / 8n z > 
d iP) = —JT = TTr TTTTTr ^ ~F P ex P 



III. SCALAR GLUEBALL DECAYS 



Since the coupling of the charmonium state to the instanton proceeds via an intermediate 
gluonic system with the quantum numbers of scalar and pseudoscalar glueballs it is natural 
to first consider direct instanton contributions to glueball decays. This problem is of course 
important in its own right. Experimental glueball searches have to rely on identifying 
glueballs from their decay products. The successful identification of a glueball requires 
theoretical calculations of glueball mixing and decay properties. In the following we compute 
the direct instanton contribution to the decay of the scalar ++ glueball state into irir, KK, 
fjf] and rjrj'. 

Since the initial state is parity even only the G 2 term in equ. (4) contributes. The relevant 
effective interaction is given by 



-I+A 



= jdz Jd (p) 

{ 



dp 



G 



X 



[(uu)(dd)(ss) + (u'j 5 u)(d'j 5 d)(ss) + (u'j 5 u)(dd)(s'j 5 s) + (uu)(d'j 5 d)(s'j 5 s)] 
(ut a u)(dt a d)(ss) + [ut a ^u)(dt a ^d)(ss) + {ut a ^u)(dt a d)(s^s) 
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+ 8L 



+ {ut a u){dt a 1 5 d){s 1 h) 
3 

--[(ut a a^u)(dt a a^d)(ss) + (ut a a^u){dt a a^d){ss) 
+ (ut a a fiul 5 u)(dt a a^d)(s 1 5 s) + (ut a a^ v u)(dt a a^d)(s^ s)} 

-^d abc [(ut a a^u)(dt b a^d)(st c s) + {ut a a^u){dt b a^d){st c s) 
+ (ut a a llul 5 u)(dt b a llu d)(st c 1 5 s) + (ut a a lxl/ u)(dt b a^ 5 d)(st c 7 5 s)\ 
+ (2 cyclic permutations u <-> d <-> s) 



-d 



abc 
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(ut a u)(dt b d)(st c s) + (ut a -? 5 u)(dt b ^d)(st c s) + (ut a ^u)(dt b d)(s^t c s) 



,5+c. 



+ {ut a u){dt b 1 5 d){st c 1 5 s) 
9 



32 



if 



abc 



(ut a a liu u)(dt a ul d)(st c a^ ll s) + (ut a a^ v Yu)(dt ov 7 o ' d)(st c a lfJ ,s) 



+ (ut a a^ 5 u) (dt b a U7 d) (st c cr 7M 7 5 s) + (ut a a^u) {dt b a ul rfd) (st c a 7 ^ 5 s) 



(5) 



Let us start with the process ++ — > tttt. In practice we have Fierz rearranged equ. (5) 
into structures that involve the strange quark condensate ss as well as operators with the 
quantum numbers of two pions. In order to compute the coupling of these operators to the 
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pions in the final state we have used PCAC relations 



(0| 57 ^.|A- + > = s K K . (7) 



The values of the decay constants are f n = 93 MeV, fx = 113 MeV [22]. We also use 
Q u = (uu) = -(248 MeV) 3 and Q d = Q u as well as Q s = 0.66Q U [23]. The coupling of the 
rj meson is not governed by chiral symmetry. A recent analysis of rj — rj' mixing and the 
chiral anomaly gives [24] 

(0\u l5 u\ri) = (0| d^d\r,) = -*(358 MeV) 2 = K*, (8) 

(0|«7 5 «|t7') = (0\d l5 d\ V ') = -i(320 MeV) 2 = K*,, 

(0|s7 5 s|^) = ^(435 MeV) 2 = K>, 

(0\s l5 s\r]'} = — z(481 MeV) 2 = K s ql . 

Finally, we need the coupling of the glueball state to the gluonic current. This quantity has 
been estimated using QCD spectral sum rules [25, 26] and the instanton model [27]. We use 

(0 ++ |(? 2 G' 2 |0) = A = 15GeV 3 . (9) 

We can now compute the matrix element for ++ — > 7r + 7r~. The interaction vertex is 

Cf + \ = Jdz Jjd (p)j^- [ (^fj (^V) 3 x -^G^ss^d^u). (10) 

The integral over the position of the instanton leads to a momentum conserving delta func- 
tion, while the vacuum dominance approximation allows us to write the amplitude in terms 
of the coupling constants introduced above. We find 

(0 ++ ( 9 )K + (p + )7r-(p-)> = (2ir) t S , (q-p + -p-)^-\ Q,Kl (11) 



dp M^ma«vt (i 2 ) 



where 

A = J p^"^> N 2_ x y a 2 j v 3 ' 
The instanton density d (p) is known accurately only in the limit of small p. For large 
p higher loop corrections and non-perturbative effects are important. The only source of 
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information in this regime is lattice QCD [28-31]. A very rough caricature of the lattice 
results is provided by the parameterization 



— = o^(p-pc), 



(13) 



p" 2 

with n ~ lfm -4 and p c ~ 0.33 fm. This parameterization gives A = (379 MeV) -9 . Another 
way to compute A is to regularize the integral over the instanton size by replacing d(p) with 
d(p) exp(— ap 2 ). The parameter a can be adjusted in order to reproduce the size distribution 
measured on the lattice. We notice, however, that whereas the instanton density scales as 
p b ~ 5 ~ p 4 , the decay amplitude scales as p b+8 ~ p 17 . This implies that the results are very 
sensitive to the density of large instantons. We note that when we study the decay of a 
small-size bound state the integral over p should be regularized by the overlap with the 
bound state wave function. We will come back to this problem in section IV below. 

We begin by studying ratios of decay rates. These ratios are not sensitive to the instanton 
size distribution. The decay rate ++ — > 7r + 7r~ is given by 



r ++„ 



1 \frn, 



o++ 



4 m 2 



16tt 



m, 



o++ 



A 

16tt 



(14) 



The decay amplitude for the process ++ — > 7r 7r is equal to the ++ — > it + it~ amplitude as 
required by isospin symmetry. Taking into account the indistinguishability of the two 7r we 
get the total nn width 

3 \jml ++ 



r ++- 



4 m 2 



32tt 



m 



o++ 



(15) 



In a similar fashion we get the decay widths for the KK, rjr), rjr]' and rj'i]' channels 

o 1 i r 



'in, 



o++ 



Am 2 K - 



16tt 



m, 



o++ 



A I 2 
^oQuK K 



To++- 



o++ 



Am 2 



>7rn 



32tt 
1 



rn 



o++ 



167T 

A 



16tt 



\ K« 2{Q S K« + (Q u + Q d )K s v ) 



(16) 
(17) 



m 2 ++ - {m r] + m v ,) 2 ][ml ++ 



m v '¥ 



16tt 



m 



x 



A 



16tt 



\ (2Q s K«K q v , + (Q u + Q d ){KlK s v , + K'K%)] 



r ++_ 



1 



o++ 



4m 2 , 



32tt 



m 



o++ 



A 



16tt 



\oK%2(Q a K% + (Q u + Q d )K% 



(18) 
(19) 



Here, KK refers to the sum of the K + K~ and K K final states. We note that in the chiral 
limit the instanton vertices responsible for im and KK decays are identical up to quark 
interchange. As a consequence, the ratio of the decay rates T ++_^ 7T7T /T ++^ K ^ is given by 
the phase space factor multiplied by the ratio of the coupling constants 



r ++. 



3 2 K 4 



r 0++ ^ 4 QIKi \ml ++ -Aml 



rn 



o++ 



4?71 

-§- = (0.193 ±0.115) 



rn 



o++ 



4 m 2 



(20) 



The main uncertainty in this estimate comes from the value of m s , which is not very accu- 
rately known. We have used m s = (140 ± 20) MeV. The ratio of nn to r/rj decay rates is not 
affected by this uncertainty, 



0.69 
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rn 



o++ 



Ami 



rn, 



o++ 



Am 2 



(21) 



In Fig.l we show the decay rates as functions of the glueball mass. We have used Aq C d = 
300 MeV and adjusted the parameter a to give the average instanton size p = 0.29 fm. We 
observe that for glueball masses m ++ > 1 GeV the KK phase space suppression quickly 
disappears and the total decay rate is dominated by the KK final state. We also note that 
for m ++ > 1.5 GeV the f]r] rate dominates over the hit rate. 

In deriving the effective instanton vertex equ. (10) we have taken all quarks to be massless. 
While this is a good approximation for the up and down quarks, this it is not necessarily 
the case for the strange quark. The m s ^ contribution to the effective interaction for ++ 
decay is given by 



P o{p) n! 



3^4 



1 TT 6 p 



¥" 3 



m s p(a s G 2 ) x 



(22) 



uu)(dd) + (w 7 5 w)(d7 5 d) + £ (ut a u)(dt a d) + (u-fh a u)(d^t a d) + 



,5j.a. 



{ua^t a u) (da^ed) 



lua. 



7 5 t a u)(d<j^H a d) 



5*a r 



There is no m s ^ contribution to the KK channel. The m s ^ correction to the other 
decay channels is 

3 \fm 



o++ 



4 m 2 



0++^7T7T — 



0++-^ — 



32tt 



m, 



o++ 



1 



o++ 



Am 2 



32tt 



m 



o++ 



167T 
1 

16tt 



X K 2 n (AQ s - 2Bm s 



(23) 



A 2[{AQ S - 2Bm s )(K*) 2 + A{Q U + Q d )K'K«)] 
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resonance 


full width T(MeV) 


Mass (MeV) 


decay channels 


/o(1370) 


200-500 


1200-1500 


pp dominant 
tttt, KK, rjrj seen 








V r KR = 0.19 ± 0.07 

-L 7T7T 


/o(1500) 


109 ±7 


1507 ± 5 


^ = 0.095 ± 0.026 

-L 7T7T 

^ = 0.18 ± 0.03 

J- irir 


/o(1710) 


125 ± 10 


1713 ±6 


$™ - 0.39 ±0.14 

1 KK 

r„„ = .4 8 ± o.l5 

1 KK 



TABLE I: Masses, decay widths, and decay channels for scalar-isoscalar mesons with masses in the 
(1 — 2) GeV range. The data were taken from [22]. 
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o++ 
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n + m r] ,) 2 }[ml ++ - {m v - m^) 2 } 
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m: 



o++ 



x 



167T 



A [2 (AQ S - 2Bm s )K«K«, + A(Q U + Q d )(K*K% + K^Kfi] 



where 



m 2 0++ - 4m 2 , 



32tt 



m: 



o++ 



16n 



X 2[(AQ S - 2Bm s )(K q rj/ ) 2 + A(Q U + Q d )K%K% 



7r 3 p 4< 
a 2 



A N 2 



(24) 



The decay rates with the m s ^ correction to the instanton vertex taken into account are 
plotted in Fig. 2. We observe that effects due to the finite strange quark are not negligible. 
We find that the 7T7T , f]f]', and ri'i]' channels are enhanced whereas the r]7] channel is reduced. 
For a typical glueball mass m ++ = (1.5 — 1.7) GeV the ratio r = B(tctt)/B(KK) changes 
from r ~ 0.25 in the case m s = to r ~ 0.55 for m s ^ 0. In Fig. 3 we show the dependence of 
the decay rates on the average instanton size p. We observe that using the phenomenological 
value p = 0.3 fm gives a total width r ++ ~ 100 MeV. We note, however, that the decay 
rates are very sensitive to the value of p. As a consequence, we cannot reliably predict the 
total decay rate. On the other hand, the ratio of the decay widths for different final states 
does not depend on p and provides a sensitive test for the importance of direct instanton 
effects. 

In Tab. I we show the masses and decay widths of scalar-isoscalar mesons in the (1-2) GeV 
mass range. These states are presumably mixtures of mesons and glueballs. This means 
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that our results cannot be directly compared to experiment without taking into account 
mixing effects. It will be interesting to study this problem in the context of the instanton 
model, but such a study is beyond the scope of this paper. It is nevertheless intriguing that 
the /o(1710) decays mostly into KK. Indeed, a number of authors have suggested that the 
/o(1710) has a large glueball admixture [32-35]. 



IV. ETA CHARM DECAYS 

The 7] c is a pseudoscalar J PC = 0~ + charmonium bound state with a mass m Vc = (2979 ± 
1.8) MeV. The total decay width of the r/ c is r^ c = (16 ± 3) MeV. In perturbation theory 
the total width is given by 

r M , = «( 1 + ^). W 

Here, i/j(0) is the 1 S ground state wave function at the origin. Using m c = 1.25 GeV and 
«s(^c) = 0.25 we get \ip{0)\ ~ 0.19GeV 3,/2 , which is consistent with the expectation from 
phenomenological potential models. Exclusive decays cannot be reliably computed in pertur- 
bative QCD. As discussed in the introduction Bjorken pointed out that r\ c decays into three 
pseudoscalar Goldstone bosons suggest that instanton effects are important [10]. The rele- 
vant decay channels and branching ratios are B(KKn) = (5.5±1.7)%, B(rjirir) = (4.9±1.8)% 
and B(t]'itit) = (4.1 ± 1.7%). These three branching ratios are anomalously large for a single 
exclusive channel, especially given the small multiplicity. The total decay rate into these 
three channels is (14.5 ±5.2)% which is still a small fraction of the total width. This implies 
that the assumption that the three-Goldstone bosons channels are instanton dominated is 
consistent with our expectation that the total width is given by perturbation theory. For 
comparison, the next most important decay channels are _B(2(7r + 7r~)) = (1.2 ± 0.4)% and 
B(pp) = (2.6 ± 0.9)%. These channels do not receive direct instanton contributions. 

The calculation proceeds along the same lines as the glueball decay calculation. Since 
the r] c is a pseudoscalar only the GG term in equ. (4) contributes. The relevant interaction 
is 

{[(u>y 5 u)(dd)(ss) + (uu)(d>y 5 d)(ss) + (uu)(dd)(s>y 5 s) + (u>y 5 u)(d>y 5 d)(s>y 5 s)] 

10 



3 
+ 8 



{uf-i s u){dt a d)(ss) + (ut"u)(dt a i s d)(3s) + (ut a u) (dfd)(tfs) 



a. ,5 , 



+ (ut a Yu)(dt a Yd)(s 1 5 s) 

3 - 
--[(ut a a^ 5 u)(dt a a^d)(ss) + (ut a a^u)(dt a a^ 5 d)(ss) 

+ (ut a a^u)(dt a a^d)(s^s) + (ut a a^u)(dt a a^d)(s^ s)} 

-^d abc [(ut a a^ 5 u)(dt b a, u d)(st c s) + (ut a a^u)(dt b a, ul 5 d)(st c s) 

+ (ut a a IMU u)(dt b a^d)(st c 1 5 s) + (^V^7 5 M )(JtV^ 7 5 rf)(st c 7 5 s)] 

+ (2 cyclic permutations u <-> d <-> s) 
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40' 



(Mt a 7 5 M)((itV)(st C s) + (Mt a w)(cr 7 5 d)(st C s) + (Mt a M)(rfrrf)(s7 5 t C s) 



+ (Mt a 7 5 M)(rftV^)(st c 7 5 s) 
9 



-—if abc 
32 J 



(ut a a tlv ') ' u)(dt a vl d){st c a llx s) + (ut a a /iu u)(dt £7„ 7 1 < d)(st c a llx s) 



(26) 



+ (u^a^u) (dt b a ul d) (st c a lfJi ^s) + (ut a a^ 5 u) (dt b a vl ^ b d) {st c a llx ^ b s) 

The strategy is the same as in the glueball case. We Fierz-rearrange the lagrangian (26) 
and apply the vacuum dominance and PCAC approximations. The coupling of the r\ c bound 
state to the instanton involves the matrix element 



\ c = ( Vc \g 2 GG\0). 



(27) 



We can get an estimate of this matrix element using a simple two-state mixing scheme for 
the rjc and pseudoscalar glueball. We write 



\rj c ) = cos(9)\cc)+sm(9)\gg), 
|(T+) = -sm(0)\cc)+cos(6)\gg). 



(28) 
(29) 



The matrix element f Vc = (0|2m c c7 5 c|?7 c ) ~ 2.8 GeV 3 is related to the charmonium wave 
function at the origin. The coupling of the topological charge density to the pseudoscalar 
glueball was estimated using QCD spectral sum rules, A -+ = (0\g 2 GG\0~ + ) ~ 22.5 GeV 3 
[26]. Using the two-state mixing scheme the two "off-diagonal" matrix elements / -+ = 
(0|2m c c7 5 c|CT + ) and \ Vc = (Ol^GGI^e) are given in terms of one mixing angle 9. We can 
estimate this mixing angle by computing the charm content of the pseudoscalar glueball 
using the heavy quark expansion. Using [36] 



C7 c = 



' a s GG + O ( — Q 



8nm 



(30) 
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we get / -+ ~ 0.14 GeV 3 and a mixing angle 9 ~ 3°. This mixing angle corresponds to 

A„ c ~ 1.12GeV 3 . (31) 

The uncertainty in this estimate is hard to assess. Below we will discuss a perturbative 
estimate of the instanton coupling to r\ c . In order to check the phenomenological consistency 
of the estimate equ. (31) we have computed the r/ c contribution to the {g 2 GG(0)g 2 GG(x)) 
correlation function. The results are shown in Fig. 4. The contribution of the pseudoscalar 
glueball is determined by the coupling constant A -+ introduced above. The couplings of the 
i], rf and ?7(1440) resonances can be extracted from the decays J/ip — > 777 [37]. We observe 
that the r] c contribution is strongly suppressed, as one would expect. We also show the 
i] c and CT + glueball contributions to the (c7 5 c(0)c7 5 c(x)) correlation function. We observe 
that even with the small mixing matrix elements obtained from equs. (28-30) the glueball 
contribution starts to dominate the rj c correlator for x > 1 fm. 

We now proceed to the calculation of the exclusive decay rates. There are four final states 
that contribute to the KKir channel, r] c -> K + K~ir°, K°K°ir°, K+K°ir- and K~K Q ir+. 
Using isospin symmetry it is sufficient to calculate only one of the amplitudes. Fierz rear- 
ranging equ. (26) we get the interaction responsible for the r\ c — > K + K~ir° 

r3„4\ ,A N 3 



C? + T*° =JdzJ poip)^^ (^j (Vp 3 )' \{a s GG){s^u){u^s){d^d). (32) 
The decay rate is given by 



T K + K - n o = J (phase space) x \M\' 



1 12 
-.A\ c K n K 2 K 



x (0.111 MeV), (33) 



.I671V2 

with A given in equ. (12). Isospin symmetry implies that the other KKir decay rates are 
given by 

^K+K-n° = ^RORO-k = (~/^) ^K f >K-n+ = ^K+K®k-- (34) 



The total KKir decay rate is 
In a similar fashion we obtain 



r - ?x 

L ryn-K 2 



I 2 

' A\ c K^K 2 K 



16ttV2 



x (0.111 MeV). (35) 



x (0.135 MeV), (36) 
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r^'7T7r 



KKr) 



KKr)' 



r)t)r) 



3 

7; x 
2 


■ 1 

.167T 


2 x 


1 

J. 

167T 


2 x 


1 

167T 


1 

- X 
6 


" 3! 

16tt 



x (0.0893 MeV) , 



A\ c K^K 2 K x (0.0788 MeV) , 



x (0.0423 MeV), 
x (0.0698 MeV) . 



(37) 
(38) 
(39) 

(40) 



Here, the first factor is the product of the isospin and final state symmetrization factors. 
The second factor is the amplitude and the third factor is the phase-space integral. 

In Fig. 5 we show the dependence of the decay rates on the average instanton size. We 
observe that the experimental KKir rate is reproduced for p = 0.29 fm. This number is 
consistent with the phenomenological instanton size. However, given the strong dependence 
on the average instanton size it is clear that we cannot reliably predict the decay rate. On 
the other hand, the following ratios are independent of the average instanton size 

2 



^KKn 

r 

- 1 - r/mr 


= 4 x 


r 


( KS r, 






r KKr) 
^KKn 


1 

= - X 

3 


^ KKr) 


_ ( Kq n 


^ KKr)' 




r 

^ KKir 


1 

" 36 X 



K 2 

K 



x 



x 



0.135 \ 
0.0893/ 



0.111 
0.135 

= 1.01, 



4.23 ± 1.27, 



V2KI 



'0.0788\ 

y 1 - "J 



x 



0.0788\ 
0.0423/ 



0.111 
2.91, 



0.141 ±0.042, 



3\y/2(K«) 2 K, 



K n K 2 K 



x 



0.0698 
0.111 



0.011 ±0.003, 



(41) 
(42) 
(43) 
(44) 
(45) 



where we have only quoted the error due to the uncertainty in m s . These numbers should 
be compared to the experimental results 



KKtt 



r)7V7V 



r)7rn 



= 1.1 ±0.5 



exp 



r»/7r7r 



= 1.2 ±0.6. 



(46) 
(47) 



exp 



We note that the ratio B(t]itit)/ B(r]'iriT) is compatible with our results while the ratio 
B(KKir)/B(r]7r7r) is not. This implies that either there are contributions other than instan- 
tons, or that the PCAC estimate of the ratio of coupling constants is not reliable, or that 
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the experimental result is not reliable. The branching ratios for r/mi and r/'mr come from 
MARK II/III experiments [38, 39]. We observe that our results for B(KKt])/B(KKtt) and 
B(KKrf')/B(KKn) are consistent with the experimental bounds. 

Another possibility is that there is a significant contribution from a scalar resonance that 
decays into nir. Indeed, instantons couple strongly to the a(600) resonance, and this state is 
not resolved in the experiments. We have therefore studied the direct instanton contribution 
to the decay i] c — > arj. After Fierz rearrangement we get the effective vertex 



C av = J dA (a s GG) - (urfu)(dd)(ss) + (uu)(drfd)(ss) + (uu)(dd)(srf s) 



I 



dB m s (a s GG) - (wy 5 u)(dd) + (uu)(dj 5 d) 
2 



(48) 



where the integrals A and B are defined in equ. (12,24). The only new matrix element we 
need is f a = (a\uu + dd\0) ~ (500 MeV) 2 [40]. We get 

1 



m vc - ( m ° + m v) 2 ][ m vc - ( m * - m v) 2 ] 



16tt 



UKUQs - 2Bm,)K« + AK*Q d ] 



(49) 



Compared to the direct decay r\ c — > rjim the r\ c — > r/cr channel is suppressed by a factor 
~ (2-n 2 / m 2 nc ) ■ (Q q f a /Kl) 2 ~ 1/100. Here, the first factor is due to the difference between 
two and three-body phase space and the second factor is the ratio of matrix elements. We 
conclude that the direct production of a a resonance from the instanton does not give a 
significant contribution to i] c — > r)(r]')7nr. This leaves the possibility that the 7T7t channel is 
enhanced by final state interactions. 

Finally, we present a perturbative estimate of the coupling of the i] c to the instanton. 
We follow the method used by Anselmino and Forte in order to estimate the instanton 
contribution to i] c — > pp [41]. The idea is that the charmonium state annihilates into two 
gluons which are absorbed by the instanton. The Feynman diagram for the process is shown 
in Fig. 6. The amplitude is given by 

d A k\ f d^ki 



A 



cc^I 



(27T) 



v(p 2 ) 



(27T) 



7m 



2 i> x - & - m c ,u 2 



«(pi) A^ d (k 2 )A b l j d (k 1 ), 



(50) 



where u(p) and v(p) are free particle charm quark spinors and A^ (k) is the Fourier trans- 
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form of the instanton gauge potential 

Air 2 fi a k v / 1 \ 

A^\k) = -*—^rm, m = 4 (1 - -K 2 {kp){k P f) . (si) 

The amplitude for the charmonium state to couple to an instanton is obtained by folding 
equ. (50) with the r) c wave function ip(p). In the non-relativistic limit the amplitude only 
depends on the wave function at the origin. 

The perturbative estimate of the transition rate is easily incorporated into the results 
obtained above by replacing the product AX Vc in equs. (33-40) according to 

AK - / po( P ) (^V) 3 (4.) ^|^(0)|/, c (p) x ffop-, (52) 

T ( , _ [i4 h fc 2 $(fc)$(fc - 2 Pc ) 

WXP) -Jdk ^ _ _ + ^ . (53 ) 

Here, p c = (m c , 0) ~ (M Vc /2,0) is the momentum of the charm quark in the charmonium 
rest frame. We note that because of the non-perturbative nature of the instanton field higher 
order corrections to equ. (52) are only suppressed by g 2 {rn~ v ) / g 2 (p)- 

The integral I Vc cannot be calculated analytically. We use the parameterization 

vr 2 A) p 4 log(l + l/(m c p)) 



with 



T ( ) u < &v 1 J. ^LLL 

" AP> - 1 + B (m c p) 4 log(l + l/(m c p))' 



(54) 



which incorporates the correct asymptotic behavior. We find that A = 0.213 and B = 0.124 
provides a good representation of the integral. In Fig. 7 we show the results for the t] c decay 
rates as a function of the average instanton size. We observe that the results are similar to 
the results obtained from the phenomenological estimate equ. (31). The effective coupling 
(A\ Vc ) differs from the estimate equ. (31) by about a factor of 3. The experimental KKn 
rate is reproduced for p = 0.31 fm. 



V. CHI CHARM DECAYS 



Another interesting consistency check on our results is provided by the study of instanton 
induced decays of the Xc into pairs of Goldstone bosons. The Xc is a scalar charmonium 
bound state with mass m Xc = 3415 MeV and width T Xc = 14.9 MeV. In a potential model 
the Xc corresponds to the 3 P state. In perturbation theory the total decay rate is dominated 
by cc — > 2g. The main exclusive decay channels are Xc 2(/T + 7r~) and Xc vr + 7r~ K + K~ 
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with branching ratios (2.4±0.6)% and (1.8±0.6)%, respectively. It would be very interesting 
to know whether these final states are dominated by scalar resonances. We will concentrate 
on final states containing two pseudoscalar mesons. There are two channels with significant 
branching ratios, Xc ~~ * vr + 7r~ and Xc —* K + K~ with branching ratios (5.0 ± 0.7) ■ 10~ 3 and 
(5.9 ±0.9) • 10~ 3 . 

The calculation of these two decay rates proceeds along the same lines as the calculation 
of the ++ glueball decays. The only new ingredient is the Xc coupling to the gluon field 
strength G 2 . We observe that the total Xc decay rate implies that (0\2m c cc\x c ) = 3.1 GeV 3 ~ 
(0\2m c ci^f 5 c\r] c ) . This suggests that a rough estimate of the Xc coupling to G 2 is given by 

A Xc = (rfG 2 |0> ~ X Vc = 1.12 GeV 3 . (55) 

Using this result we can obtain the Xc decay rates by rescaling the scalar glueball decay 
rates equ. (23-24) according to 

2 



^Xc^ml,m2 — ro++-»ml,m2 X f ^ ^ 



(56) 

m 0++ -+m Xc 



where ml,m2 labels the two-meson final state. In Fig. 8 we show the dependence of the Xc 
decay rates on the average instanton size p. We observe that the experimental 7r + 7r~ decay 
rate is reproduced for p = 0.29 fm. In Fig. 9 we plot the ratio of decay rates for 7r + 7r~ and 
K + K~ . Again, the experimental value is reproduced for p ~ 0.3 fm. 

Finally, we can also estimate the cc coupling to the instanton using the perturbative 
method introduced in section IV. In the case of the Xc we use 

1 ^ a . 1 rrrvtf* f rf °(p) , /4 2 3 \VK) 



^ - i^^'(0) / ^ (rV) >< 'M (57) 



4tt Xc 2v^F V X y ' J P 5 ' V3 ' ) g 2 (p) 

4tt Xc 2^ V X J P 5 V3 11 g 2 {p) 



where R'(0) — 0.39 GeV 5 ^ 2 is the derivative of the 3 P wave function at the origin and I Xc 
is the loop integral 

r (^ - f^u Hk)H\*Pc-k\) 15(k- Pc ) 2 + 3m 2 c + 4k 2 

In Fig. 10 we compare the perturbative result with the phenomenological estimate. Again, 
the results are comparable. The experimental 7r + 7r~ rate is reproduced for p = 0.29 fm. 
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VI. SUMMARY 



In summary we have studied the instanton contribution to the decay of a number of 
"gluon rich" states in the (1.5-3.5) GeV range, the scalar glueball, the r/ c and the Xc- In the 
case of charmonium instanton induced decays are probably a small part of the total decay 
rate, but the final states are very distinctive. In the case of the scalar glueball classical fields 
play an important role in determining the structure of the bound state and instantons may 
well dominate the total decay rate. 

We have assumed that the gluonic system is small and that the instanton contribution to 
the decay can be described in terms of an effective local interaction. The meson coupling to 
the local operator was determined using PCAC. Using this method we find that the scalar 
glueball decay is dominated by the KK final state for glueball masses m ++ > 1 GeV. In the 
physically interesting mass range 1.5 GeV < m ++ < 1.75 GeV the branching ratios satisfy 
B(ryq) : B(nn) : B(KK) = 1 : (3.3 ± 0.3) : (5.5 ± 0.5). 

Our main focus in this work are i] c decays into three pseudoscalar Goldstone bosons. 
We find that the experimental decay rate T(r] c — > KKir) can be reproduced for an average 
instanton size p = 0.31, consistent with phenomenological determinations and lattice results. 
This in itself is quite remarkable, since the phenomenolgical determination is based on 
properties of the QCD vacuum. 

The ratio of decay rates B(r)'7nr) : B(i]7i7i) : B(KKtt) = 1:1: (4.2 ± 1.3) is insensitive 
to the average instanton size. While the ratio B(r)'irir) : B(r)inr) = 1 : 1 is consistent with 
experiment, the ratio B(rjirir) : B(KK) = 1 : (4.2±1.3) is at best marginally consistent with 
the experimental value 1.1 ± 0.5. We have also studied Xc decays into two pseudoscalars. 
We find that the absolute decay rates can be reproduced for p = 0.29 fm. Instantons are 
compatible with the measured ratio B(K + K~) : B(tt + tt~) = 1.2 

There are many questions that remain to be answered. On the experimental side it would 
be useful if additional data for the channels r\ c — > rj'inr, r\ixix were collected. One important 
question is whether (nn) resonances are important. It should also be possible to identify 
the smaller decay channels r\ c — > KKr], KKr/'. In addition to that, it is interesting to study 
the distribution of the final state mesons in all three-meson channels. Instantons predict 
that the production mechanism is completely isotropic and that the final state mesons are 
distributed according to three-body phase space. 
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In addition to that, there are a number of important theoretical issues that remain to be 
resolved. In the limit in which the scalar glueball is light the decay ++ — > itit(KK) can be 
studied using effective lagrangians based on broken scale invariance [42-44]. Our calculation 
based on direct instanton effects is valid in the opposite limit. Nevertheless, the instanton 
liquid model respects Ward identities based on broken scale invariance [16] and one should 
be able to recover the low energy theorem. In the case ++ — > mr(KK) one should also be 
able to study the validity of the PCAC approximation in more detail. This could be done, 
for example, using numerical simulations of the instanton liquid. Finally we need to address 
the question how to properly compute the overlap of the initial cc system with the instanton. 
This, of course, is a more general problem that also affects calculations of electroweak baryon 
number violation in high energy pp collisions [45, 46] and QCD multi-particle production in 
hadronic collisions [47]. 
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Glueball decay rates vs glueball mass, corresponding to p average =0.29 fm and m s =0 MeV 
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FIG. 1: Scalar glueball decay rates plotted a function of the mass of the scalar glueball. The rates 
shown in this figure were computed from the instanton vertex in the chiral limit. The average 
instanton size was taken to be p = 0.29 fm. 
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Glueball decay rates vs glueball mass, corresponding to p aV erage = - 29 fm > m s =140 MeV 
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FIG. 2: Same as Fig. 1 but with m s / corrections in the instanton vertex taken into account. 
The results shown in this figure correspond to m s = 140 MeV. 
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Glueball decay rates vs P average for m ++=1 600 MeV and m s =140 MeV 
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FIG. 3: Dependence of glueball decay rates on the average instanton size. The results shown in 
this figure correspond to the instanton vertex with m s / terms included. The strange quark 
mass was taken to be m s = 140 MeV. 
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FIG. 4: Resonance contributions to the pseudoscalar glueball correlation function 
(g 2 GG(0)g 2 GG(x)) and the charmonium correlator (0750(0)0750(0;)) . Both correlation functions 
are normalized to free field behavior. In the case of the gluonic correlation function we show the 
glueball contribution compared to the 77, rj', r/(1440) and rj c contribution. For the charmonium 
correlation function we show the r\ c and glueball contribution. 
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r| c main decay rates vs average instanton size 
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FIG. 5: Decay widths r) c — > KKir and rf c — > r/iTir as a function of the average instanton size p. The 
short dashed line shows the experimental KKir width. 




FIG. 6: The Feynman diagram corresponding to the perturbative treatment of charmonium decay. 
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r| c perturbative and phenomenological decay rates vs average instanton size 
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FIG. 7: Decay rates T{r] c — > KKtt) and T(r] c — > ry7T7r) as a function of the average instanton size p. 
We show both the results from a phenomenological and a perturbative estimate of the cc coupling 
to the instanton. 
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Chi decay rates in instanton-specific channels vs p 
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FIG. 8: Decay widths \c — ► K + K , ir + ir and ??ry as a function of the average instanton size p. 
The short dashed line shows the experimental K + K~ width. 
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Ratio of decay rates for % vs p average 
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FIG. 9: Ratio B(\c — ► tt + tt~)/ 'B(xc K + K~) of decay rates as a function of the average 
instanton size. The dashed line shows the experimental value 0.84. We also show the experimental 
uncertainty, as well as the uncertainty in the instanton prediction due to the the value of the 
strange quark mass. 



28 



X decay rates in instanton-specific channels vs p average for perturbative and phenomenological methods 
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FIG. 10: Decay rates T(x c — ► vr + 7r _ ) and T(x c —> K + K~) as a function of the average instanton 
size p. We show both the results from a phenomenological and a perturbative estimate of the cc 
coupling to the instanton. 
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